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Abstract. We consider the semi-infinite optimization problem: 
/* := min {/(x) : 5 (x, y) < 0, Vy G Y x }, 

xEEX 

where f,g are polynomials and X C M™ as well as Y x C W, x e X, are 
compact basic semi-algebraic sets. To approximate /* we proceed in two steps. 
First, we use the "joint+marginal" approach of the author [9] to approximate 
from above the function x h-> "^(x) = sup{g(x, y) : y £ Y x } by a polynomial 
&d > ^5 °f degree at most 2d, with the strong property that ^ converges to 
$ for the Li-norm, as d — > oo (and in particular, almost uniformly for some 
subsequence (di), £ S N). Therefore, to approximate /* one may wish to solve 
the polynomial optimization problem /9 = min xe x{/(x) : < 0} via 

a (by now standard) hierarchy of semidefinite relaxations, and for increasing 
values of d. In practice d is fixed, small, and one relaxes the constraint &d < 
to $(j(x) < e with e > 0, allowing to change e dynamically. As d increases, the 
limit of the optimal value /j is bounded above by /* + e. 



1. Introduction 

Consider the semi-infinite optimization problem: 

(1.1) P: f* :=min{/(x):p(x,y) < 0, Vy G Y x }, 

xex 

where X C M n , Y x C W for every x G X, and some functions / : M n — > R, 
j:i"xK p :4 E. 

Problem P is called a semi-infinite optimization problem because of the infin- 
itely many constraints g(x, y) < for all y G Y x (for each fixed x G X). It has 
many applications and particularly in robust control. 

In full generality P is a very hard problem and most methods aiming at com- 
puting (or at least approximating) /* use discretization to overcome the difficult 
semi-infinite constraint g(x, y) < for all y G Y x . Namely, in typical approaches 
where Y x = Y for all x G X (i.e. no dependence on x), the set Y C l p is dis- 
cretized on a finite grid and if the resulting nonlinear programming problems 
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are solved to global optimality, then convergence to a global optimum of the 
semi-infinite problem occurs as the grid size vanishes (see e.g. the discussion 
and the many references in [IQ]). Alternatively, in pU] the authors provide lower 
bounds on /* by discretizing Y and upper bounds via convex relaxations of the 
inner problem max ye Y{5'(x, y)} < 0. In [11] the authors also use a discretiza- 
tion scheme of Y but now combined with a hierarchy of sum of squares convex 
relaxations for solving to global optimality. 

Contribution. We restrict ourselves to problem P where : 

• f,g are polynomials, and 

• XcM" and Y x C MP, x G X, are compact basic semi-algebraic sets. 

For instance many problems of robust control can be put in this framework; see 
e.g. their description in [J]. Then in this context we provide a numerical scheme 
whose novelty with respect to previous works is to avoid discretization of the set 
Y x . Instead we use the "joint +marginal" methodology for parametric polynomial 
optimization developed by the author in [H] , to provide a sequence of polynomials 
C R[x] (with degree 2d, d G N) that approximate from above the function 
$(x) := max y {g(x, y) : y G Y x }, and with the strong property that if d — > oo 
then <3>d — > $ in the Li-norm. (In particular, $^ — > $ almost uniformly on X 
for some subsequence (de), £ G N.) Then, ideally, one could solve the nested 
sequence of polynomial optimization problems: 

(1.2) P d . / d * = min{/(x) : $ d (x) < }, d=l,2,... 

For fixed d, one may approximate (and often solve exactly) (ll.2p by solving 
a hierarchy of semidefinite relaxations, as defined in [6]. However, as the size 
0(d n ) of these semidefinite relaxations increases very fast with d, in practice one 
rather let d be fixed, small, and relax the constraint 3>d( x ) < to $d(x) < e 
for some scalar e > that one may adjust dynamically during the algorithm. 
As d increases, the resulting optimal value f\ is bounded above by /* + e. The 
approach is illustrated on a sample of small problems taken from the literature. 

2. Notation, definitions and preliminary results 

Let R[x] (resp. K.[x,y]) denote the ring of real polynomials in the variables 
x = (x 1 ,...,x n ) (resp. x and y = (y 1 , . . . ,y p )), whereas £[x] (resp. E[x,y]) 
denote its subset of sums of squares. 

Let K[y]fe C R[y] denote the vector space of real polynomials of degree at most 
k. For every a G N n the notation x a stands for the monomial x" 1 ■ ■ ■ x^™ and 
for every d G N, let := {a G : £ < d} with cardinal s(d) = ( n + d ). 
Similarly N p d := {0 G N p : £\ ^ < d} with cardinal ( p+ p d ). A polynomial 
/ 6 R[x] is written 

x->/(x) = f" x °> 
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and / can be identified with its vector of coefficients f = (f a ) in the canonical 
basis. For a real symmetric matrix A the notation A y stands for A is positive 
semidefinite. 

A real sequence z = (z a ), a G N n , has a representing measure if there exists 
some finite Borel measure \x on W 1 such that 

z a = I x a d/i(x), VaG N n . 

Given a real sequence z = (z a ) define the linear functional L z : R[x] — > K. by: 
/(=^/ Q x a ) ^L z (/) = £/ a z Q , /£%]. 

Moment matrix. The moment matrix associated with a sequence z = (z a ), 
a G N n , is the real symmetric matrix M^(z) with rows and columns indexed 
by N^, and whose entry (a, /3) is just z a +p, for every a, f3 G N^. If z has a 
representing measure /x then M d (z) y because 

(f,M d (z)f) = J fdfjL > 0, Vf G R'W. 

Localizing matrix. With z as above and g G M[x] (with g(x) = X) 7 #7 x7 )> ^e 
localizing matrix associated with z and g is the real symmetric matrix NLd(gz) 
with rows and columns indexed by N^, and whose entry (a, (3) is just 9f z a+p+y, 
for every a, (3 G N^. If z has a representing measure \x whose support is contained 
in the set {x : g(x) > 0} then M^gz) y because 

(f,M^z)f) = J fgdfi >0, Vf Gl s(rf) . 

Definition 2.1 (Archimedean property). A set of polynomials qj G M[x], j = 
0, . . . ,p (with q Q = 1), satisfy the Archimedean property if the quadratic polyno- 
mial x i—)- M — ||x|| 2 can be written in the form: 

p 

^- l|x|| 2 = ^^-(x)g,(x), 

3=0 

for some sums of squares polynomials (aj) C S[x]. 

Of course the Archimedean property implies that the set D := {x G R n : 
9j( x ) > 0, j — 1, . . . ,p} is compact. For instance, it holds whenever the level 
set {x : (7fc(x) > 0} is compact for some k G {1, . . . or if the q/s are affinc 
and D is compact (hence a polytope). On the other hand, if D is compact then 
M — ||x| 2 > for all x G D and some M sufficiently large. So if one adds the 
redundant quadratic constraint x h-> g p+ i(x) = M — ||x|| 2 > in the definition of 
D then the Archimedean property holds. Hence it is not a restrictive assumption. 
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Let D := {x G K™ : 5 3 -(x) > 0, j — 1, . . . , p}, and given a polynomial /i G 
consider the hierarchy of semidefmite programs: 



(2.1; 



Pi = min L z (h) 

s.t. M e (z),M e _ V] {q,z) h0, j = l,...,p, 

where z = (z a ), a G N% e , and u,- = |~(deggj)/2], j — 1, . . . ,p. 

Theorem 2.2 (j6j [8]). Lei a family of polynomials (qj) C K.[x] satisfy the 
Archimedean property. Then as i — > oo, pe f h* = min x {/i(x) : x G D}. 
Moreover, if z* is an optimal solution of $2. 1\) and 

(2.2) rankM^z*) = rankM f _„(z*) (=: r) 

(where v = maxjVj) then pi = h* and one may extract r global minimizers 
x* G D ; k = l,...,r. 

The size (resp. the number of variables) of the semidefmite program (12.11) 
grows as ( n ^ £ ) (resp. as (™^ 2£ )) and so becomes rapidly prohibitive, especially 
in view of the present status of available semidefmite solvers. Therefore, and 
even though practice reveals that convergence is fast and often finite, so far, the 
above methodology is limited to small to medium size problems (typically, and 
depending on the degree of the polynomials appearing in the data, problems with 
up to n G [10,20] variables). However, for larger size problems with sparsity in 
the data and/or symmetries, adhoc and tractable versions of (12.11) exist. See for 
instance the sparse version of (12. lj) proposed in [12] . and whose convergence was 
proved in [7] when the sparsity pattern satifies the so-called running intersection 
property. In [12] this technique was shown to be successful on a sample of non 
convex problems with up to 1000 variables. 

3. Main result 

Let B C M™ be a simple set like a box or an ellipsoid. Let p s G M[x], s = 
1, . . . , sx, and hj G R[x, y], j — 1, . . . , m, be given polynomials and let X C M n 
be the basic semi-algebraic set 

X := {x G R n : p s (x) > 0, s = l,...,sx}. 

Next, for every xeR", let Y x C M p be the basic semi-algebraic set described by: 

(3.1) Y x = {y eW : ^(x,y) > 0, j = l,...,m}, 
and with B D X, let K C R n x W be the set 

(3.2) K := {(x, y) G R n+P : x G B; ^(x,y)>0, j = l,...,m}. 
Observe that problem P in (II. ip is equivalent to: 

(3.3) P: /* = min{/(x) : $(x) < 0} 

(3.4) where $(x) = max{g(x, y) : y G Y x }, x G B. 
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Lemma 3.1. Let K C M. n+P in A3.2\) be compact and assume that for every 
x G B C M. n , the set Y x defined in / Iff, ij) is nonempty. Then $ is upper semicon- 
tinuous (u.s.c.) on B. Moreover, if there is some compact set YcK f such that 
Y x = Y for every x G B ; then $ is continuous on B. 

Proof. Let xo G B be fixed, arbitrary, and let {xk)k£N C B be a sequence that 
converges to x and such that 

limsup $(x) = lim $(xfc). 

x— s-x A;— >oo 

As K is compact then so is Y x for every x G B. Therefore, as Y x 7^ for 
all x G B and g is continuous, there exists an optimal solution yj* G Y Xfe for 
every k. By compactness there exist a subsequence (hi) and y* G M p such that 
(x^,y*J (x ,y*) G K, as £ -> 00. Hence 

limsup $(x) = lim $(xfe) 



/[■- 



= lim #(x fe ,y£) = lim #(x fe ,,y£ ) 

fc— >-oo t— >oo 

= #(x ,y*) < $(x ), 

which proves that $ is u.s.c. at x . As x G B was arbitrary, $ is u.s.c. on B. 

Next, assume that there is some compact set Y C MP such that Y x = Y for 
every x G B. Let x G B be fixed arbitrary with $(x ) = g(x ,yo) for some 
yo G Y. Let (x n ) C B, n G N, be a sequence such that x n — > x as n — )■ 00, 
and $(x ) > liminf $(x) = lim $(x n ). Again, let y* G Y be such that $(x„) = 

x -x, 11 >:x. 

g(x n ,y*), n G N. By compactness, consider an arbitrary converging subsequence 
(ni) C N, i.e., such that (x nf ,y*J -> (x ,y*) G K as £ -> 00, for some y* G Y. 
Suppose that $(x ) (= ^(xo,yS)) > #(xo,y*), say $(x ) > g(x ,y*) + 5 for some 
5 > 0. By continuity of g, g(x nv y * ) < ^(xo, y*) + 5/2 for every £ > £\ (for some 
ti). But again, by continuity, ^(x^yjj) — p(x ,yo)| < 5/3 whenever £ > £ 2 (for 
some £ 2 )- And so we obtain the contradiction 

$(x n J > <7(x n „y*) > $(xo)-5/3 
$(x n J = #(x^,y*J < $(xo)-5/2, 

whenever £ > max[li,^]- Therefore, g(x ,yo) = g(x ,y*) and so, 

#(x ,yo) = $(xo) = #(x ,y*) = lim $(x n J = liminf $(x) < $(x ), 

£— >oo x— »xo 

which combined with $ being u.s.c, yields that $ is continuous at xo. □ 

We next explain how to 

• approximate the function x 1— > $(x) on B by a polynomial, and 

• evaluate (or at least approximate) $(x) for some given x G B, to check 
whether $(x) < 0. 

Indeed, these are the two main ingredients of the algorithm that we present later. 

5 



3.1. Certificate of $(x) < 0. For every x G X fixed, let g XJ h* G M[y] be the 
polynomials y i— >■ g x (y) = g(x,y) and y i— >■ hj(y) := fy(x,y), j = 1, . . . ,m, and 
consider the hierarchy of semidefinite programs: 



(3.5) Q,(x) : 



p t (x) = max L z (a x ) 

it. M e (z),M e - Vj (hJz) to, j = l, 



where z = (zp), (3 G N^, and Vj = |~(deg h*)/2~\ , j — 1, . . . ,m. Obviously one has 
P£( x ) > ^( x ) f° r every £, and 

Corollary 3.2. Lei x G X and assume that the polynomials (/i*) C M[y] satisfy 
the Archimedean property. Then: 

(a) v4s £ — > oo, pt{x) 4 ^( x ) = max{#(x, y) : y G Y x }. In particular, if 
P^( x ) < /or some i, then $(x) < 0. 

(b) Moreover, if z* an optimal solution of \3. 51) £na£ satisfies 

rankM £ (z*) = rankM £ _ t) (z*) (=: r), 

(where v := max,,- u,,), inen p^(x) = $(x) and inere are r global maximizers 
y(k) G Y x; fc = l,...,r. 

Corollary 13.21 is a direct consequence of Theorem 12.21 



3.2. Approximating the function $. Recall that B D X is a simple set like 
e.g., a simplex, a box or an ellipsoid and let p be the finite Borel probability 
measure uniformly distributed on B. Therefore, the vector 7 = (ja), a G N n , of 
moments of p, i.e., 

7a := I x a dp(x), a G N n , 
Jb 

can be computed easily. For instance, in the sequel we assume that B = [—1, l] n = 
{x : #j(x) > 0, i = l,...n} with $i G R[x, y] being the polynomial (x, y) 1— > 
6»i(x,y) := 1 - ac?, % = 1, . . .,n. 

Observe that the function $ is defined in H3.4[) via a parametric polynomial 
optimization problem (with x being the parameter vector). Therefore, following 
[9], let Tj = |~(deg hj)/2] , j — 1, . . . , m, and consider the hierarchy of semidefinite 
relaxations indexed by d G N: 



(3.6) 



max L z (g) 

s.t. M d (z),M d _ r .(fyz) ^0, j = l,...,m 
M dr _ l (9 i z) h0, z = l,...,n 



L z (x Q ) = 7q , a G N 



n 

2d- 



where the sequence z is now indexed in N^" p , i.e., z = (z a p), (a, (3) G 

Writing g = 1, the dual of the semidefmite program ( 13. 6 p reads 

(3.7) 

P*d = mir i / g(x) dfi(x) 

m n 

s.t. g(x) -^(x,y) = J^o- i (x,y)/i i (x,y)+^V<(x,y)0i(x,y) 

g G R[x] 2(i , c.pV'i e S[x,y] 
deg (jj /ij < 2d, j = 0, . . . , m. 
deg^i^j < 2d, z = 

It turns out that any optimal solution of the semidefmite program (13.71) permits 
to approximate $ in a strong sense. 

Theorem 3.3 (|9J). Let K C M. n+P in VS.ty be compact. Assume that the polyno- 
mials hj, 0i G ]R[x, y] satisfy the Archimedean property and assume that for every 
x G B ; the set Y x defined in A3.1\) is nonempty. Let $d G M[x]2d be an optimal 
solution of \3. 7| ). Then : 

(a) $d > $ and as d — > oo ; 

(3.8) / ($ d (x)-$(x))a>(x) = I | $ d (x) - $(x) | d//(x) ->• 0, 

i/ioi is, $d — t- $ /or £/ie Li(B, [i)-norm^. 

(b) T/iere zs a subsequence (de), £ G N, suc/i £/iai $^ — >■ $, ^-almost uniformly^ 
in B ; as £ — >■ oo. 

The proof of (a) can be found in j9], whereas (b) follows from (a) and [U The- 
orem 2.5.3]. 



3.3. An algorithm. The idea behind the algorithm is to approximate P in (II. ip 
with the polynomial optimization problem: (P^): 

(3.9) P^: /J = min{/(x) : $ d (x) < e }, d=l,2,... 

xex 

with d G N, e > fixed, and as in Theorem I3.3[ for every d = 1, . . .. 

Obviously, for e = one has /° > /* for all d because by definition $d > $ 
for every d G N. However, it may happen that P° has no solution. Next, if x* 
is an optimal solution of P and $(x*) < 0, it may also happen that $d(x*) > 
if d is not large enough. This is why one needs to relax the constraint $ < 
to $d < e for some e > 0. However, in view of Theorem I3.3[ one expects that 
fd ~ /* provided that d and e are sufficiently large and small, respectively. And 
indeed: 

1 Li(B,/i) is the Banach space of /i-integrable functions on B, with norm ||/|| = f B \ f\dfj,. 
2 If one fixes e > arbitrary then there is some A £ 6(B) such that fJ-(A) < e and <&d e & 
uniformly on B \ A, as £ — > oo. 
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Theorem 3.4. Assume that X is the closure of an open set. Let e > be fixed, 
arbitrary and with f e d be as in / Iff. 9\) , let G X be any optimal solution of |5J|) 
(including the case where e = 0), and let 

Id '■= min {/l : Z = l,...,d} = /(x| (d) ) f or some £(d) G {1, . . . , d}. 

(a) If e > i/iere exists d t &N such that for every d > d t , /(x^) < f* + e. 

(b) // there is an optimal solution x* G X o/ (COP suc/i i/iat $(x*) < ; then 
there exists do G N such that for every d > d , f* < /(x^) < /* + e. 

Proof, (a) With e > fixed, arbitrary, let x* G X be such that $(x*) < and 

/(x*) < /* + e/2. We may assume that x* is not on the boundary of X. Let 

0\ := {x G int X : $(x) < e/2} which is an open set because $ is u.s.c. (by 

Lemma [37TT) . and so p{0\) > 0. Next, as / is continuous, there exists po > 

such that /</* + £ whenever x 6 O f 2 := {x 6 int X : ||x — x*|| < po}- Observe 

that p := p(0\ H O e 2 ) > because 0\ n O e 2 is an open set (with x* G O, 1 n O 2 ). 

Next, by Theorem 13.3( b). there is a subsequence (o^), £ G N, such that $^ —> $, 

//-almost uniformly on B. Hence, there is some Borel set A e C B, and integer 

£ e G N, such that p(A e ) < p/2 and sup |$(x) - $^(x)| < e/2 for all £ > 4- In 

xex\A e 

particular, as p{A e ) < p/2 < //(O^flO 2 ), the set A e := (0^n0 2 )\A e has positive 
//-measure. Therefore, /(x) < /* + e and $d £ (x) < e whenever £ > £ t and x G A e , 
which in turn implies /J f < /* + e, and consequently, /J = /(xj/^) < /* + e, the 
desired result. 

(b) Let e' := -$(x*), and let 0\, := {x G int X : $(x) < -e'/2} which is a 
nonempty open set because it contains x* and $ is u.s.c. Let O 2 be as O 2 in 
the proof of (a), but now with x*, = x* G X. Both 0\, and O 2 are open and 
nonempty because they contain x*. The rest of the proof is like for the proof of 
(a), but noticing that now for every x G A e / one has $^ ( x ) < — e '/2 + e'/2 = 0, 
and so x is feasible for (13.91) with e = 0. Next, by feasiblity /(x) > /* since the 
resulting feasible set in (13.91) is smaller than that of (11.21) because $^ > $, for 
all d. And so /* < /(x) < /* + e whenever x G A e , and £ > £ e , from which (b) 
follows. □ 

Theorem 13.41 provides a rationale behind the algorithm that we present below. 
In solving (13.91) with d sufficiently large and small e (or even e = 0), /j would 
provide a good approximation of /*. But in principle, computing the global 
optimum /j is still a difficult problem. However, P e d is a polynomial optimization 
problem. Therefore, by Theorem 12.21 if the polynomials (p s ) C R[x] that define 
X satisfy the Archimedean property (see Definition 12. ip we can approximate /J 
from below, as closely as desired, by a monotone sequence (/jj, t G N, obtained 
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by solving the hierarchy of semidefinite relaxations ( 12. ip . which here read: 

( fa = min L z (f) 

( 3 - 10 ) I s.t. M,(z),M,_ d (e-$ d z) ^0 

{ M t _ ts (p s z) ^ 0, s = l,...,sx, 

where t s = |~(degp s )/2] , s = 1, . . . , sx. 

Corollary 3.5. Assume that the polynomials (p s ) C R[x] satisfy the Archimedean 
property. Then f e dt "f /j as t — >• oo. Moreover, ifz* is an optimal solution of A3.10\) 
and 

(3.11) rankM^z*) = rankM t _t (z*) (=: r) 

(where to := max[d, max s [i s ]]j i/ien /J t = /j anc? one may extract r global mini- 
mizers x^(&) G X, = 1, . . . , r. 77mt is, /or every k = 1, . . . , r, f(x^(k)) = f d 
and $ d (x.* d (k)) < e. 

However, given a minimizer x. d 6 X, if on the one hand $^(x^) < e, on the other 
hand it may not satisfy $(x^) < 0. (Recall that checking whether $(x^) < can 
be done via solving the hierarchy of relaxations Q^(x) in (13. 5p with x := x^.) If 
this happens then one solves again ( 13. 10[) for a smaller value of e, etc., until one 
obtains some x^ G X with $(x^) < 0. 

Finally, and as already mentioned, if d is relatively large, the size of semidefinite 
relaxations ( I3.10P to compute f dt becomes too large for practical implementation 
(as one must have t > d). So in practice one let d be fixed at a small value, 
typically the smallest possible value of d, i.e., 1 (<3>d is quadratic) or 2 ($^ is 
quartic)), and one updates e as indicated above. So the resulting algorithm 
reads: 

Algorithm. 

Input: £, d, k* G N, e > (e.g. e := 1(H), d G N, x := *, /(*) = +oo. 

Output: /(xj) with x^ G X and $(xj) < 0. 

Step 1: Set = 1 and e{k) = 1. 

Step 2: While fc < fc*, solve P^ fc) in (EHJ) -»■ x* G X. 

Step 3: Solve Q £ (x*) in ->• M x fc)- 

If -e < p^ x fe) < set x* := x* and STOP. 

If P£(x* k ) < -e Q then: 

• if /(x) > /(x£) then set x := x£. If k = k* then xjj := x£. 

• set e(fc + 1) := 2e(fc), k :— k + 1 and go to Step 2. 
If p £ (x*) > then: 

• If k < k* set e(k + 1) := e(k)/2, k := k + 1 and go to Step 2. 

• If A; = k* then set set x^ = x. 

Observe that in Step 2 of the above algorithm, one assumes that by solving 
P e } k) one obtains x£ G X. 
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3.4. Numerical experiments. We have taken Examples 2, 7, 9, K, M, N, all 
from Bhattacharjee et al. [2J Appendix A] and whose data are polynomials, ex- 
cept for problem L. For the latter problem, the non-polynomial function x i— > 
min[0, (xi — x 2 )} is semi-algebraic and can be generated by introducing an addi- 
tional variable £3, with the polynomial constraints: 

x\ = (xi - x 2 ) 2 ] x 3 > 0. 

Indeed, 2min[0, (x\ — x 2 )} = X\ — x 2 — x 3 . 

Although these examples are quite small, they are still non trivial (and even 
difficult) to solve, and we wanted to test the above methodology with small 
relaxation order d. In fact we have even considered the smallest possible d, i.e., 
d — 1 is quadratic). Results in Table 1 are quite good since by using the 
semidefinite relaxation of minimal order "cf' one obtains an optimal value /J 
quite close to /*, at the price of updating e several times. 

Next, for Problem L, if we now increase d to d — 2, we improve the opti- 
mal value which becomes /J = 0.3849 with e = 2.2. However, for Problem M, 
increasing d does not improve the optimal value. 
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best known value 


fd 


final value of e 


problem 2 


0.194 


0.198 


1.895 


problem 7 


1.0 


1.41 


5 


problem 9 


-12.0 


-14.47* 





problem K 


-3.0 


-3.0 


3.037 


problem L 


0.3431 


0.435 


2.295 


problem M 


1.0 


2.25 


2.592 


problem N 


0.0 


10" 8 






Table 1. Examples of [2> Table 6.1] with minimal d 



4. Conclusion 

We have presented an algorithm for semi-infinite (global) polynomial optimiza- 
tion whose novelty with respect to previous works is to not rely on a discretization 
scheme. Instead, it uses a polynomial approximation $^ of the function $, ob- 
tained by solving some semidefinite relaxation attached to the "joint +marginal" 
approach developed in [9] for parametric optimization, which guarantees (strong) 
convergence <3>d — » $ in Li-norm. Then for fixed d, one has to solve a polynomial 
optimization problem, which can be done by solving an appropriate hierarchy 
of semidefinite relaxations. Of course, as already mentioned and especially in 
view of the present status of semidefinite solvers, so far the present methodology 
is limited to small to medium size problems, unless sparsity in the data and/or 
symmetries are taken into account appropriately, as described in e.g. [7J IL2] . 
Preliminary results on non trivial (but small size) examples are encouraging. 
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